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Summary. The high sensitivity of JWST will open a new window on the end of the cosmo- 
logical dark ages. Small stellar clusters, with a stellar mass of several x 10''Mp,, and low-mass 
black holes (BHs), with a mass of several x IO^Mq should be directly detectable out to red- 
shift z = 10, and individual supernovae (SNe) and gamma ray burst GRB afterglows are bright 
enough to be visible beyond this redshift. Dense primordial gas, in the process of collaps- 
ing from large scales to form protogalaxies, may also be possible to image through diffuse 
recombination line emission, possibly even before stars or BHs are formed. In this article, I 
discuss the key physical processes that are expected to have determined the sizes of the first 
star-clusters and black holes, and the prospect of studying these objects by direct detections 
with JWST and with other instruments. The direct light emitted by the very first stellar clusters 
and intermediate-mass black holes at z > 10 will likely fall below JWST?, detection thresh- 
old. However, JWST could reveal a decline at the faint-end of the high-redshift luminosity 
function, and thereby shed light on radiative and other feedback effects that operate at these 
early epochs. JWST^iW also have the sensitivity to detect individual SNe from beyond z = 10. 
In a dedicated survey lasting for several weeks, thousands of SNe could be detected at z > 6, 
with a redshift distribution extending to the formation of the very first stars at z > 15. Using 
these SNe as tracers may be the only method to map out the earliest stages of the cosmic star- 
formation history. Finally, we point out that studying the earliest objects at high redshift will 
also offer a new window on the primordial power spectrum, on ~ 100 times smaller scales 
than probed by current large-scale structure data. 



1 Introduction 

The formation of the first astrophysical objects and the subsequent epoch of reionization are 
at the frontiers of research in Astronomy. Recent years have seen significant progress both 
in our theoretical understanding and in observational probes of this transition epoch in the 
early universe. Through a combination of methods, including measurements of the cosmic 
microwave background (CMB) anisotropies, culminating in the precise determination of the 
temperature and polarization power spectra from the Wilkinson Microwave Anisotropy Probe 
(WMAPj experiment 1 16 159 , 160 , 134 42 , 96 1, Hubble diagrams of distant Supernovae |739l 
II38I , and various probes of large scale structure (see, e.g., references in |5] |I72||T56l ) the key 
cosmological parameters have been determined to high accuracy. 
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Because of the emergence of a concordance (ACDM) cosmology, we can securely pre- 
dict the collapse redshifts of the first non-linear dark matter condensations: 2 — 3(T peaks of 
the primordial density field on mass scales of 10^^^ M,;), corresponding to the cosmological 
Jeans mass, collapse at redshifts z= 15 — 20. Of course, in addition to the cosmological pa- 
rameters describing the average background universe (ACDM), one also needs a description 
of the seed fluctuations, to make such predictions. The nearly scale-invariant nature of the ini- 
tial fluctuation power spectrum has also been empirically confirmed by the above-mentioned 
large-scale structure data. However, the predictions for the earliest halos rely on extrapolat- 
ing the primordial power spectrum to a wavenumber of A: ~ 10 /iMpc^', a scale that is 2-3 
orders of magnitude smaller than the smallest scale directly probed by current data (e.g. 11781 
and references therein). It is possible that the small-scale power differs substantially from the 
extrapolated value. For example, in the case of warm dark matter (WDM), the power can be 
reduced by many orders of magnitude on the relevant scales, and, to zeroth order, the first gen- 
eration of halos predicted in A CDM would not exist 1 6 1 . Recent simulations indicate that the 
formation of the first halos would indeed be delayed I i l32l|193l , and the details of the collapse 
dynamics modified, possibly ultimately affecting the properties of the first stars [61 ]. Within 
the ACDM paradigm, however, robust predictions can be made, using three-dimensional sim- 
ulations I195II1621 . rather than semi-analytical methods f 137"154V Such predictions are now 
limited mainly by the ~ 5% uncertainty in the normalization of the primordial power spec- 
trum, tTg/,-i (e.g. |T56|), which translates essentially to a ~ 5% uncertainty in the collapse 
redshift of the first structures. 

On the observational side, the recent measurement of the optical depth to electron scatter- 
ing (Tp = 0.09 ± 0.03) by WMAP suggests that the first sources of light significantly ionized 
the intergalactic medium (IGM) at redshift z ~ 11 ±3 1 134 160 1. The Sloan Digital Sky Sur- 
vey (SDSS) has uncovered a handful of bright quasars at redshifts as high as z = 6.41 (see 
1 141 1461 l46l 1471 l48l |49J or the recent review by 1501 ). The luminosity of these sources rivals 
those of the most powerful quasars at the peak of their activity at z — 2.5. A sizable popula- 
tion of z ~ 6 galaxies have also been identified in the past several years (see, e.g., |44 | for a 
recent comprehensive review), many of them in broad-band photometry in deep, wide-area 
Hubble Space Telescope fields Ii35lll 12I|19II164||165I . Narrow-band searches have also been 
successful in discovering numerous z > 5 Lyman-a emitting galaxies. Starting with the initial 
success of 11821 . 10m class telescopes have now identified about two dozen such objects out 
to z ~ 6.7 (see 1 1 1 1] for a compilation of datasets; and the more recent Subaru survey results 
by, e.g., 1 169 93J). In summary, the above observational data leave little doubt that the first 
galaxies and quasars were indeed well in place by redshift z — 6, but not in significant numbers 
prior to redshift z — 15. 

These developments have been complemented by progress in theoretical work focusing on 
the cooling and collapse of gas in primordial low-mass halos (see reviews by 1 7 1 , 2 1 and ref- 
erences therein). A broad picture has emerged, identifying several processes that are important 
in the formation of the first stars and black holes (BHs) out of the baryonic gas in the earliest 
halos. The key ingredient in this picture is the abundance of H2 molecules that form via gas- 
phase reactions in the early universe (this beautifully simple fact was recognized as early as in 
1967 by 1 147 1). The first objects form out of gas that cools via H2 molecules, and condenses 
at the centers of virialized dark matter "minihalos" with virial temperatures of Tvir ~ 200K 
1821 11711 . Detailed numerical simulations have shown convergence toward a gas tempera- 
ture T ~ 300 K and density n ~ 10^ cm^^, dictated by the thermodynamic properties of H2 
fTll20l[T9? 194 1, which allows the collapse of a single clump of mass 10^-10^ Mq at the cen- 
ter of the halo. The soft UV and X-ray radiation emitted by the stars and black holes formed in 
the first handful of such clumps provide prompt and significant global feedback on the chem- 
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leal and thermal state of the IGM. A relatively feeble early background radiation field can al- 
ready have strong effects on gas cooling and H2 chemistry, and can produce significant global 
impact, affecting the formation of the bulk of the first astrophysical structures, and therefore 
the reionization history of the IGM |[80l|72l[l06l[l07l[l29l[l^|97][5iE^ The ob- 

servational data summarized above suggests that such feedback processes indeed shaped the 
reionization history, and hints that star-formation in early minihalos was suppressed i73i . 



2 Summary of Relevant Physics at High Redshifts 

In the context of ACDM cosmologies, it is natural to identify the first dark matter halos as the 
hosts of the first sources of light. A simplified picture for the emergence of the first sources 
of light, and consequent reionization is as follows: the gas in dark matter halos cools and 
turns into ionizing sources (stars and/or black holes), which produce UV (and possibly X- 
ray) radiation, and drive expanding ionized regions into the IGM. The volume filling factor of 
ionized regions then, at least initially, roughly tracks the global formation rate of dark halos. 
Eventually, the ionized regions percolate (when the filling factor Fmi reaches unity), and the 
remaining neutral hydrogen in the IGM is rapidly cleared away as the ionizing background 
builds up. However, a soft UV radiation background at photon energies of < 13.6eV, as well 
as possibly a soft X-ray background at > IkeV, will build up before reionization, since the 
early IGM is optically thin at these energies. 

We now briefly describe the various effects that should determine the evolution of the first 
sources and of reionization at high redshifts. For in-depth discussion, we refer the reader to 
extended reviews (e.g. ||9l|), and to reviews focusing on the roles of H2 molecules in reioniza- 
tion (2), on the effect of reionization on CMB anisotropies |75 |, and on progress in the last 
three years in studying reionization |71 , 54|. 

In the discussions that follow, it will be useful to distinguish three different types of dark 
matter halos, which can be roughly divided into three different ranges of virial temperatures, 
as follows: 

300K < Tvir < 10"* K (Type II - susceptible to Ha-feedback) 

10^ K < Tyif < 2 X 10^ K (Type la - susceptible to photo-heating feedback) 

Tyir > 2 X 10^ K (Type lb - gas can fall in, even in the face of photo-heating) 

We will hereafter refer to these three different types of halos as Type II, Type la, and Type 
lb halos. The motivation in distinguishing Type II and Type I halos is based on H2 molecular 
vs atomic H cooling, whereas types "a" and "b" reflect the ability of halos to allow infall 
and cooling of photoionized gas. While the actual values of the critical temperatures, and the 
sharpness of the transition from one category to the next is uncertain (see discussion below), as 
we argue below, star formation in each type of halo is likely governed by different physics, and 
each type of halo likely plays a different role in the reionization history. In short. Type II halos 
can host ionizing sources only in the neutral regions of the IGM, and only if H2 molecules 
are present; Type la halos can form new ionizing sources only in the neutral IGM regions, but 
can do so by atomic cooling, irrespective of the H2 abundance, and Type lb halos can form 
ionizing sources regardless of the H2 abundance, and whether they are in the ionized or neutral 
phase of the IGM. It is useful to next provide a summary of the various important physical 
effects that govern star-formation and feedback in these halos. 
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Table 1. Collapse redshifts and masses of 1, 2 and 3cT dark matter halos with virial temper- 
atures of Tyii- = 100 and lO'^K, reflecting the minimum gas temperatures required for cooling 
by H2 and neutral atomic H, respectively. 



T'virCK) V ZcoU Mh,io(M,3) 



100 


1 


7 


2x 10^ 




2 


16 


5x 10^ 




3 


26 


3x lO'* 


104 


1 


3.5 


4x 10** 




2 


9 


1 X 10** 




3 


15 


6x 10^ 



H2 Molecule Formation. In ACDM cosmologies, structure formation is bottom-up: the 
earliest nonlinear dark matter halos form at low masses. Gas contracts together with the 
dark matter only in dark halos above the cosmological Jeans mass, Mj ~ 10'*. However, 
this gas can only cool and contract to high densities in somewhat more massive halos, with 
M > Mh2 = 1O^M0[(1 +z)/ll]"-'/^ (i.e. Type II halos), and only provided that there is a suf- 
ficient abundance of H2 molecules, with a relative number fraction at least nH2/«H ~ 10^^ 
1821 11711 . Because the typical collapse redshift of halos is a strong function of their size, 
the abundance of H2 molecules is potentially the most important parameter in determining 
the onset of reionization. For example, 2(7 halos with virial temperatures of lOOK appear at 
z = 16, while 2a halos with virial temperatures of lO^K (Type la halos, in which gas cooling 
is enabled by atomic hydrogen lines) appear only at z = 9. Table [T] summarizes the collapse 
redshifts and masses of dark halos at these two virial temperatures in the concordance cos- 
mology. As a result, the presence or absence of H2 in Type II halos makes a factor of ~2 
difference in the redshift for the onset of structure formation, and thus potentially effects the 
reionization redshift and the electron scattering opacity t by a similar factor. In the absence of 
any feedback processes, the gas collecting in Type II halos is expected to be able to form the 
requisite amount of H2 182111711 . However, both internal and external feedback processes can 
alter the typical H2 abundance (see discussion below). 

The Nature of the Light Sources in the First Halos. A significant uncertainty is the na- 
ture of the ionizing sources turning on inside halos collapsing at the highest redshifts. Three 
dimensional simulations using adaptive mesh refinement (AMR; 1 1 1) and smooth particle hy- 
drodynamics (SPH; |20, 194|) techniques have followed the contraction of gas in Type II halos 
at high redshifts to high densities. These works have shown convergence toward a tempera- 
ture/density regime of T ~ 200 K, n ^ lO'^cm^"', dictated by the critical density at which the 
excited states of H2 reach equilibrium population levels. These results have suggested that the 
first Type II halos can only form unusually massive stars, with masses of at least 100 Mq, 
and only from a small fraction (< 0.01) of the available gas. Such massive stars are effec- 
tive producers of ionizing radiation, making reionization easier to achieve. In addition, it is 
expected that the first massive stars, forming out of metal-free gas, have unusually hard spec- 
tra I176II22II148I , which is important for possibly ionizing helium, in addition to hydrogen. 
An alternative possibility is that a similar fraction of the gas in the first Type II halos forms 
massive black holes L76.,110.,144J . In fact, when massive, non-rotating metal-free stars end 
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their life, they are expected to leave behind stellar-mass seed BHs, unless their mass is in the 
range of 140-260 M0 [84|. These early black holes can then accrete gas (possibly only after 
a delay - if the progenitor star clears the host halo of gas, the BH will have to await until the 
host halo merges with another, gas-rich halo), acting as "miniquasars". The miniquasars will 
produce a hard spectrum extending to the soft X-rays, which could be important in catalyzing 
H2 formation globally (see discussion below). 

Efficiencies and the Transition from Metal Free to "Normal" Stars. Another fundamental 
question of interest is the efficiency at which the first sources inject ionizing photons into the 
IGM. This can be parameterized by the product — Ny /*/esc) where f^ — j (^^b^halo I ^m) 
is the fraction of baryons in the halo that turns into stars; Ny is the mean number of ionizing 
photons produced by an atom cycled through stars, averaged over the initial mass function 
(IMF) of the stars; and /esc is the fraction of these ionizing photons that escapes into the 
IGM. It is difficult to estimate these quantities at high-redshifts from first principles, but the 
discussion below can serve as a useful guide. 

Although the majority of the baryonic mass in the local universe has been turned into 
stars L56i , the global star formation efficiency at high redshifts was likely lower. To explain a 
universal carbon enrichment of the IGM to a level of 10^^ — 10^"^ Zq, the required efficiency, 
averaged over all halos at z > 4, is /* = 2 — 20% |I761. However, the numerical simulations 
mentioned in § 2.2 above suggest that the fraction of gas turned into massive stars in Type II 
halos is/, < 1%. 

The escape fraction of ionizing radiation in local starburst galaxies is of order ~ 10%. The 
higher characteristic densities at higher redshifts could decrease this value I41II187I , although 
there are empirical indications that the escape fraction in z ~ 3 galaxies may instead be higher 
[166|, at least in some galaxies 1 153 1. Radiation can also more readily ionize the local gas, and 
escape from the small Type II halos which have relatively low total hydrogen column densities 
(< 10'^ cm^^), effectively with /esc = 1 in the smallest halos 1 183 1. 

The ionizing photon yield per proton for a normal Salpeter IMF is Ny ~ 4000. However, 
if the IMF consists exclusively of massive M > 2OOM0 metal-free stars, then Ny can be up 
to a factor of ~ 20 higher 1221 11481 . The transition from metal-free to a "normal" stellar 
population is thought to occur at a critical metallicity of Zcr ~ 5 x 1O"^Z0, above which 
cooling and fragmentation becomes efficient and stops the IMF from being biased toward 
massive stars 1221 . It is natural to associate this transition with that of the assembly of halos 
with virial temperatures of > lO'^K (Type la halos). Type II halos are fragile, and likely blow 
away their gas and "shut themselves off" after a single episode of (metal-free) star-formation. 
They are therefore unlikely to allow continued formation of stars with metallicities above Zerit. 
Subsequent star-formation will then occur only when the deeper potential wells of Type la 
halos are assembled and cool their gas via atomic hydrogen lines. The material that collects 
in these halos will then have already gone through a Type II halo phase and contain traces of 
metals. 

As we argue in §|3]below, there exists an alternative, equally plausible scenario. Most Type 
II halos may not have formed any stars, due to global H2 photodissociation by an early cos- 
mic soft~UV background. In this case, the first generation of metal-free stars must appear in 
Type la halos. Halos above this threshold can eject most of their self-produced metals into the 
IGM, but, in difference from Type II halos, can retain most of their gas U081 , and can have 
significant episodes of metal-free star formation. These halos will also start the process of 
reionization by driving expanding ionization fronts into the IGM. The metals that are ejected 
from Type la halos will reside in these photoionized regions of the IGM. As discussed in §[3] 
below, photoionization heating in these regions may suppress gas infall and cooling, causing 
a pause in the formation of new structures, until larger dark matter halos, with virial tempera- 
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tures of Tyir > 2 x lO^K (Type lb halos) are assembled. The material that collects in Type lb 
halos will then have already gone through a previous phase of metal-enrichment by Type la 
halos, and it is unlikely that Type lb halos can form significant numbers of metal-free stars. 



3 Feedback Effects 

Several feedback effects are likely to be important for modulating the evolution of the early 
cosmic star formation history and reionization. There can be significant internal feedback in 
or near each ionizing source, due to the presence of supernovae 1531 , or of the radiation field 
I131II1421 . on the local H2 chemistry. The net sign of these effects is difficult compute, as it 
depends on the source properties and spectra, and on the detailed density distribution internal 
and near to the sources. For practical purposes of computing the feedback that modulates the 
global star-formation or reionization, we may, however, think of any internal feedback effect 
as regulating the efficiency parameter e,, defined above. 

Since the universe is optically thin at soft UV (below 13.6eV), and soft X-ray (>lkeV) 
photon energies, radiation from the earliest Type II halos can build up global backgrounds at 
these energies, and provide prompt external, global feedback on the formation of subsequent 
structures, which are easier to follow. A very large number of studies over the past several 
years have assessed the feedback from both the LW and X-ray backgrounds quantitatively, 
including many works that employed three-dimensional simulations 1 72 32 1061 11071 11421 
|I43ll65]|97ll5llllI6l |4lFlJ. In isolation. Type II halos with virial temperatures as low as 
a few 100 K could form enough H2, via gas-phase chemistry, for efficient cooling and gas 
contraction I82II17II . However, H2 molecules are fragile, and can be dissociated by soft UV 
radiation absorbed in their Lyman-Werner (LW) bands |80 32l |72| [T4T1 . In patches of the 
IGM corresponding to fossil HII regions that have recombined (after the death of short-lived 
ionizing source, such as a massive star), the gas retains excess entropy for a Hubble time. This 
"entropy floor" can reduce gas densities in the cores of collapsing halos (analogously to the 
case of "preheating" of nearby galaxy clusters), and decrease the critical LW background flux 
that will photodissociate H2 1 129|. On the other hand, positive feedback effects, such as the 
presence of extra free electrons (beyond the residual electrons from the recombination epoch) 
from protogalactic shocks 1152! 1531 , from a previous ionization epoch II291 11681 . or from 
X-rays |80lll27i[T42l|I43l , can enhance the H2 abundance. 

The extent to which star-formation in Type II halos was quenched globally has remained 
unclear, with numerical simulations generally favoring less quenching |106, 107, 142, 143] 
1 1 861 than predicted in semi-analytical models. Likewise, simulations find a smaller, if any, 
effect from X-rays |97|. In a recent study |116|, we used three-dimensional hydrodynamic 
simulations to investigate the effects of a transient ultraviolet (UV) flux and a LW background 
on the collapse and cooling of pregalactic clouds, with masses in the range 10^ - 10^ Ma,, at 
high redshifts (z > 18). In the absence of a LW background, we found that a critical specific 
intensity of iuv ~ 0.1 (in units of 10^^'ergs s^' cm^^ Hz^' sr^') demarcates the transition 
from net negative to positive feedback for the halo population. Note that this flux is ~2 orders 
of magnitude below the level required for reionization (defined by the requirement of produc- 
ing a few photons per H atom). A weaker UV flux stimulates subsequent star formation inside 
the fossil HII regions, by enhancing the H2 molecule abundance. A stronger UV flux signifi- 
cantly delays star-formation by reducing the gas density, and increasing the cooling time at the 
centers of collapsing halos. At a fixed yuv. the sign of the feedback also depends strongly on 
the density of the gas at the time of UV illumination. In either case, once the UV flux is turned 
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olf, its impact starts to diminish after ~ 30% of the Hubble time. In the more realistic case 
when a permanent LW background is present (in addition to a short-lived ionizing neighbor), 
with 7lw <^0-01 X 10^^'ergs s^' cm^^ Hz^' sr^', strong suppression persists down to the 
lowest redshift (z = 18) in these simulations. The feedback was also found to depend strongly 
on the mass of the Type II halo, with the smaller halos strongly suppressed, but the larger 
halos (> 10^ M0) nearly immune to feedback. In recent works, 1 1 861 and 1133] found that 
the largest Type II halos, with Tvi,- > 4000K, can eventually cool their gas even in the face 
of a LW background with Ju^ > 0.1, due to the increased electron abundance and elevated 
temperature and cooling rate in the inner regions of these halos. 

A second type of important feedback is that photo-ionized regions are photo-heated to a 
temperature of > W^K, with a corresponding increase in the Jeans mass in these regions, and 
possible suppression of gas accretion onto low-mass halos (e.g. I43lll73ll67lll50l ). Reioniza- 
tion is then expected to be accompanied by a drop in the global SFR, corresponding to a sup- 
pression of star formation in small halos (i.e. those with virial temperatures below Tyir < 10** 
- 10^ K). The size of such a drop is uncertain, since the ability of halos to self-shield against 
the ionizing radiation is poorly constrained at high redshifts. Early work on this subject, in 
the context of dwarf galaxies at lower redshift 1 173|, suggested that an ionizing background 
would completely suppress star formation in "dwarf galaxy" halos with circular velocities 
^circ ^ 35 km s^', and partially suppress star-formation in halos with 35 km s^' < I'circ ^£ 
100 km s^'. However, more recent studies 1951 1361 find that at high-redshifts (z > 3), self- 
shielding and increased cooling efficiency could be strong countering effects. These calcula- 
tions, however, assume spherical symmetry, leaving open the possibility of strong feedback 
for a halo with non-isotropic gas profile, illuminated along a low-column density line of sight 
(see 11511 for a detailed treatment of three-dimensional gas dynamics in photo-heated low- 
mass halos). 

Because the earliest ionizing sources formed at the locations of the rare density peaks, 
their spatial distribution was strongly clustered. Since most feedback mechanisms operate over 
a limited length scale, their effects will depend strongly on the spatial distribution of halos 
hosting ionizing sources. Numerical simulations are a promising way to address feedback 
among clustered sources, since they capture the full, three-dimensional relationships among 
the host halos |89|. However, the dynamic range required to resolve the small minihalos, 
within a large enough cosmic volume to be representative, remains a challenge, especially in 
simulations that include radiative transfer 1881 . Semi-analytical models avoid the problems 
associated with the large dynamical range; they are also an efficient way to explore parameter 
space and serve as important sanity checks for more complicated simulations. Semi-analytical 
studies to date have included either various feedback effects (e.g. 1 771 1741 [1891 1281 1701 1581 
11881 1911 ) or the effect of source clustering on the HII bubble-size distribution (e.g. 1591 ), 
but not both. In a recent study 1981 . we have incorporated photo-ionization feedback, in a 
simplified way, into a model that partially captures the source clustering (i.e., only in the radial 
direction away from sources). Source clustering was found to increase the mean HII bubble 
size by a factor of several, and to dramatically increase the fraction of minihalos that are 
suppressed, by a factor of up to --^ 60 relative to a randomly distributed population. We argue 
that source clustering is likely to similarly boost the importance of a variety of other feedback 
mechanisms. (This enhanced suppression can also help reduce the electron scattering optical 
depth Tg, as required by the three-year data from WMAP L73j .) 
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4 How can we detect this feedback? 

There are several ways, at least in principle, to discover the presence of global feedback mech- 
anisms that modulate the early star-formation rate and reionization. Here I simply list several 
possibilities, in order to give an (admittedly crude) overview. In § (6] I will discuss in more 
detail the one most likely to be relevant to JWST - tracing the cosmic star-formation history 
with SNe, and searching for a feature in this ultra-high redshift version of the "Lilly-Madau" 
diagram, caused by the feedback. 

In general, the feedback processes discussed above can produce an extended and com- 
plex (possibly even non-monotonic) reionization history, especially at the earliest epochs. 
The earliest stages of the global reionization history (ionized fraction versus cosmic time) can 
be probed in 21cm studies (see 1571 and the contributions by Steve Fmianetto and Avi Loeb 
in these proceedings). Additionally, the measurement of polarization anisotropies by Planck 
can go beyond a constraint on the total electron scattering optical depth Tf that is measured 
by WMAP, and at least distinguish reionization histories that differ significantly 1921 1741 [T25l 
from each other. Reionization may modify the small-angle CMB anisotropies, as well, through 
the kinetic SZ effect, at a level that may be detectable in the future |145llll5iril4l|90|[T46l . 

A recent study f 1231 considered the pumping of the 63.2/lm fine-structure line of neutral 
oxygen in the high-redshift intergalactic medium (IGM), in analogy with the Wouthuysen- 
Field effect for the 21 cm line of cosmic HI. This showed that the soft UV background at 
1300A can affect the fine-structure population levels in the ground state of 01. If a signifi- 
cant fraction of the IGM volume is filled with "fossil H II regions" that have recombined and 
contain neutral OI, then this can produce a non-negligible spectral distortion in the cosmic mi- 
crowave background (CMB). A measurement of this signature can trace the global metallicity 
at the end of the dark ages, prior to the completion of cosmic reionization, and is complemen- 
tary to the cosmological 21cm studies. In addition to the mean spectrum, fluctuations in the OI 
pumping signal may be detectable, provided the background fluctuations on arcminute scales, 
around 650GHz, can be measured to the njy level 1 124|. If the IGM is polluted with metals 
at high-redshift, then CMB angular fluctuations due to density fluctuations alone (without 
pumping, just due to the geometrical effect of scattering) could also be detectable for several 
metal and molecular species (131. 

At lower redshift, the details of the later stages of reionization can be probed in more 
detail, through studying the statistics of Lyman line absorption in the spectra of quasars (e.g. 
|117||118l and references therein), gamma-ray burst afterglows, and galaxies (see, e.g., ISTl 
and references therein). In particular, the drop in the cosmic star-formation history near the 
end-stages of reionization at z ~ 6 — 7 due to photo-heating could be detected in the high- 
redshift extension of the 'Lilly - Madau' diagram |T04 , 109 1, by directly counting faint galax- 
ies | 8|. In practice, the low-mass galaxies susceptible to the reionization suppression are faint 
and may fall below JWSTs detection limit. Whether or not these galaxies will be detected 
(and in sufficient numbers so that they can reveal the effect of reionization), depends crucially 
on the feedback effects discussed above, the redshift of reionization, the size of the affected 
galaxies and their typical star-formation efficiencies, as well as the amount of dust obscura- 
tion. Alternatively, by analyzing the Lyman a absorption spectra of SDSS quasars at z ~ 6, 
(29| suggested, from the non-monotonic evolution of the mean IGM opacity, that we may 
already have detected a drop in the SFR at z ~ 6. In order to improve on this current, low 
signal-to-noise result, deep, high-resolution spectra of bright quasars would be required from 
beyond the epoch of reionization at z > 6 (this could be possible with JWST, see |79 1). A sup- 
pression of low-mass galaxies would also increase the effective clustering of the reionizing 
sources (since higher-mass halos are more strongly clustered), and increase the fluctuations 



The First Stars and Black Holes 



9 



in the Lyman a forest opacity at somewhat lower redshifts. This effect may already have been 
observed 11901 . although the modeling details still matter, and the presently measured scatter 
in opacity could perhaps still be consistent with density fluctuations alone 1 103 ]. 

In §|6] below, I will return to the issue of feedback, and discuss tracing the cosmic star- 
formation history with distant SNe. 



5 Direct Detections of the First Sources 

The most basic question to ask is whether JWST could directly detect the first stars or black 
holes. As discussed above, the very first star may have formed in isolation in a 10^ Mq dark 
matter halo at z > 15, and it will then be beyond the reach of direct imaging even by JWST. 
The critical mass for detection with JWST depends on the IMF and star-formation efficiency. 
At z = 10, assuming ~ 10% of the gas in a halo turns into stars, with a normal Salpeter IMF, 
a InJy broad-band threshold at near-IR wavelengths would allow the detection of a stellar 
cluster whose mass is a several x 10^ M0 (761 . It will help if the IMF is biased toward more 
massive, and more luminous stars. If the stars were all metal-free, on the other hand, they 
would have a lower flux than metal-enriched stars at JWSTs wavelengths, due to their high 
effective temperatures which shifts their flux to (observed) UV wavelengths H76II22III48I . 
The conclusion is that while JWST is very unlikely to detect the first individual stars directly, it 
will most likely directly measure the luminosity function of faint galaxies, extending down to 
sufficiently small sizes, corresponding to the halo masses M = 10^"'° Me. This wifl directly 
probe the feedback effects discussed above. For example, a clear turn-over in the LF at the 
luminosity corresponding to the atomic cooling threshold (or lack of it) would be compelling 
evidence for H2-feedback (or lack of it). 

Likewise, one can ask whether JWST can see the first BHs directly ? The relevant threshold 

- i.e. the lowest BH mass - will again depend on the spectrum and luminosity (in terms of, 
say, the Eddington value) of the BHs. With the average spectrum of quasars at lower redshift 
1451 , and assuming Eddington luminosity, at the InJy threshold, JWST could detect a BH 
whose mass is a several x 10^ Mq (77). If the first BHs are the remnants of massive, metal- 
free stars, then their initial masses will be below this threshold. Furthermore, it is not clear 
when such stellar-seed BHs would start shining at a significant fraction of the Eddington 
limit (since the progenitor star may clear their host halo of gas, as mentioned above). On 
the other hand, if H2-suppression prevents most Type II halos from forming stars, then the 
first stars and BHs would appear in large numbers only inside more massive dark halos, with 
virial temperatures exceeding T^^- > lO'^K. How gas cools, condenses, and fragments in such 
halos is presently not well understood. Using semi-analytical toy models, [1281 argued that 
the initial cooling by atomic H allows the gas to begin to collapse - even in the face of a 
significant LW background. If the gas remained at ~ 10'*K, the high Jeans mass ~ lO*" Mq 
in these halos would suggest that a supermassive black hole (SMBH) of a similar mass may 
form at the nucleus 1 121 [24l [1791 [l5l|T05!. Such BHs could be directly detectable by JWST 
at z ~ 10, provided they shine at > 10% of their Eddington luminosity. 

Predictions for the number counts of high redshift galaxies and quasars at JWSTs thresh- 
olds at near-infrared wavelengths have been made using simple semi-analytic models I76II77I . 
Surface densities as high as several sources per square arcminute are predicted from z > 5, 
with most of these sources at z > 10. These predictions, obtained from the simplest models, 
represent the DM halo abundance, multiplied with optimistic star/BH formation efficiencies, 
and corrected for stellar/quasar duty cycles. As such, they could be considered as upper bounds 

- feedback effects will certainly reduce the counts. If the galaxies occupy a fair fraction (~ 5%) 
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of the virial radius of their host halos, then a large fraction (> 50%) of them can potentially 
be resolved with JWSTs planned angular resolution of '--^ 0.06" 1781 171. 

In addition to broad-band searches, one may look for the earliest light-sources in emis- 
sion lines. The strongest recombination lines of H and He from 5 < z < 20 will fall in the 
near-infrared bands of JWST and could be bright enough to be detectable. Specific predic- 
tions have been made for the source counts in the Ha emission line ||127|| and for the three 
strongest Hell lines that could be powered either by BHs or metal-free stars with a hard spec- 
trum I130II1771 . The key assumption is that most of the ionizing radiation produced by the 
miniquasars is processed into such recombination lines (rather than escaping into the IGM). 
Under this optimistic assumption (which is likely violated at least in the smallest minihalos 
(j_83J), the lines are detectable for a fiducial 10^ M© miniquasar at z = 10, or for a "micro- 
galaxy" with a star-formation rate of IMpyr^' . The simultaneous detection of H and He lines 
would be especially significant. As already argued above, the hardness of the ionizing contin- 
uum from the first sources of ultraviolet radiation plays a crucial role in the reionization of the 
IGM. It would therefore be very interesting to directly measure the ionizing continuum of any 
z > 6 source. While this may be feasible at X-ray energies for exceptionally bright sources, 
the absorption by neutral gas within the source and in the intervening IGM will render the 
ionizing continuum of high redshift sources inaccessible to direct observation out to ~ IkeV. 
The comparison of Ha and Hell line strengths can be used to infer the ratio of Hell to HI 
ionizing photons, Q = ^^J^'V^kJir ^ measurement of this ratio would shed light on the nature 
of the first luminous sources, and, in particular, it could reveal if the source has a soft (stellar) 
or hard (AGN-like) spectrum. Note that this technique has already been successfully applied 
to constrain the spectra of sources in several nearby extragalactic HII regions |63|. An alter- 
native method to probe the hardness of the spectrum of the first sources would be to study the 
thickness of their ionization fronts I T?nil97ll99l ; this may be possible for bright quasars in 
Lyman line absorption | 99| or in 21cm studies 11971119^ 1991 . 

Provided the gas in the high redshift halos is enriched to near-solar levels, several molec- 
ular lines may be visible. In fact, CO has already been detected in the most distant z = 6.41 
quasar [180] , and, in fact, spatially resolved [ISIJ . The detectability of CO for high redshift 
sources in general has been considered by 11551 and by 1681 . For a star formation rate of 
> 30 Mq I yr, the CO lines are detectable at all redshifts z = 5 — 30 by the Millimeter Array 
(the redshift independent sensitivity is due to the increasing CMB temperature with redshift), 
while the Atacama Large Millimeter Array (ALMA) could reveal fainter CO emission. The 
detection of these molecular lines will provide valuable information on the stellar content and 
gas kinematics in the earliest halos. 

Finally, the baryons inside high-redshift halos with virial temperatures T > 10'*K need 
to cool radiatively, in order to condense inside the dark matter potential wells, even before 
any stars or black holes can form. The release of the gravitational binding energy, over the 
halo assembly time-scale, can result in a significant and detectable Lya flux |8fl|52|. At the 
limiting line flux ~ 10^'*^ erg s^' cm^^ asec^^ of JWST, several sufficiently massive halos, 
with velocity dispersions CT > 1 20 km s^ ' , would be visible per 4' x 4' field. The halos would 
have characteristic angular sizes of ~ 10", would be detectable in a broad-band survey out to 
z ~ 6 — 8, and would provide a direct probe of galaxies in the process of forming. They may 
be accompanied by He+ Lya emission at the 10% level 81911 . but remain undetectable at 
other wavelengths. The main challenge, if such blobs are detected, without any continuum 
source, will likely be the interpretation - as is the case for the currently detected extended 
Lyman a sources (there are currently ~three dozen extended Lya blobs known, with l/3rd of 
such objects in the largest sample 11131 consistent with being powered by cooling radiation). 
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Monte Carlo calculations of Lya radiative transfer through optically thick, spherically 
symmetric, collapsing gas clouds were presented in 1371 1381 . with the aim of identifying a 
clear diagnostic of gas infall (a similar effort in 3D is being carried out 11701 ). These represent 
simplified models of proto-galaxies that are caught in the process of their assembly. Such 
galaxies can produce Lya flux over an extended solid angle, either from a spatially extended 
Lya emissivity, or from scattering effects, or both. We presented a detailed study of the effect 
of the gas distribution and kinematics, and of the Lya emissivity profile, on the emergent 
spectrum and surface brightness distribution. The emergent Lya spectrum is typically double- 
peaked and asymmetric. The surface brightness distribution is typically flat, and the detection 
of a strong wavelength dependence of its slope (with preferential flattening at the red side of 
the line) would be a robust indication that Lya photons are being generated (rather than just 
scattered) in a spatially extended region around the galaxy, as in the case of a cooling flow. 
An alternative, robust diagnostic for scattering is provided by the polarization of an extended 
Lyman a source [39 |. Spectral polarimetry (in particular, the wavelength-dependence of the 
polarization) can differentiate between Lya scattering off infalling gas and outflowing gas. 



6 Tracing The Beginning of the Cosmic Star-Formation History 
With Supernovae 

As argued above, the very first stellar clusters may fall below the direct detection of JWST. A 
promising way (and probably the only way) to directly observe the first stars individually is 
through their explosions, either as Supernovae or gamma ray bursts (GRBs). Before we present 
the expectations for SNe in detail, we first briefly discuss an alternative, GRB afterglows. In 
the years leading up to the launch of the Swift satellite^ it has been increasingly recognized 
that distant gamma ray bursts (GRBs) offer a unique probe of the high redshift universe. In 
particular, GRBs are the brightest known electromagnetic phenomena in the universe, and 
can be detected up to very high redshifts (e.g. 1 184, 101 , 33|), well beyond the redshift z ~ 
6.5 of the most distant currently known quasars |50| and galaxies |44|. There is increasing 
evidence that GRBs are associated with the collapse of short-lived, massive stars, including 
the association of bursts with star-forming regions (e.g. [55], a contribution of supernova 
light to the optical afterglow (e.g. I17II62I ). and most directly, association with a supernova 
fT63][86l. 

As a result, the redshift distribution of bursts should follow the mean cosmic star- 
formation rate (SFR). Several studies have computed the evolution of the expected GRB rate 
under this assumption, based on empirical models of the global SFR ll74lll75lll84|[Toni33l . 
Determinations of the cosmic SFR out to redshift z ~ 5 (e.g. I26ll60ll64l ) have shown that star- 
formation is already significant at the upper end of the measured redshift range, with > 10% 
of all stars forming prior to z = 5, which would result in a significant population of GRBs at 
these redshifts. Further associating star-formation with the formation rate of non-linear dark 
matter halos, and using theoretical models based on the Press-Schechter formalism |T37 |, it 
is possible to extrapolate the SFR and obtain the GRB rates expected at still higher redshifts 
I23|[31 1201. These studies have concluded that a significant fraction (exceeding several per- 
cent) of GRBs detected at Swift's sensitivity should originate at redshifts as high as z > 10. 
The spectra of bright optical/IR afterglows of such distant GRBs can provide information both 
on the progenitor, and also reveal absorption features by neutral hydrogen in the intergalactic 
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medium (IGM), and can serve as an especially clean probe of the reionization history of the 
universe ll21l|101| [3T] [T00l[T0l . 

In summary, the main advantage of GRBs is that they are bright and can be seen essen- 
tially from arbitrary far away. The main drawback, however, is that they are rare. Even the 
most optimistic among the above-listed models predict only a handful of detectable GRBs per 
year from z > 6 - indeed, in the past two years, only one GRB has been discovered at z > 6 
(GRB050904 at z = 6.3; (94)). 

6.1 Supernovae From The First Stars 

Although SNe are not as bright as GRB afterglows, they are bright compared to galaxies at 
the very faint end of the luminosity function, and individual core-collapse SNe would still be 
visible from beyond z = 10. Furthermore, since they occur much more frequently than GRBs, 
and they remain visible for a longer time than GRB afterglows, they will likely offer a better 
(and possibly only) chance to trace out the cosmic star formation history accurately, with a 
statistically significant sample extending beyond z > 10. The above is true even if the bulk of 
early supernovae are similar to the core-collapse SNe in the local universe (note that Type la 
supernovae will not have time to form at z ^ 6). In addition, the pair-instability supernovae 
from massive, metal-free stars are expected to be much brighter than Type II supernovae from 
normal (metal-enriched) stars 1851 . 

There have been several studies of the expected early supernova rate (SNR), calibrated to 
the observed metallicity of the Lyman a forest 1 122| or the observed SFR at lower redshifts 
1341 . The expected rate of pair-instability SNe from very high-z from the first generation of 
metal-free stars was studied by 11851 . In a recent work 11191 , we constructed the expected 
detection rate of high-z SNe in SNe surveys for JWST. We also quantified the prospects of 
detecting a drop in the SN rate due to photo-heating feedback at reionization. Given that SNe 
may be our only hope to directly map out the beginning of the cosmic star-formation history, 
most of the rest of this article is devoted to discussing this possibility in detail. 

6.2 The Global Star Formation and Supernova Rates 

The global SFR density can be obtained and extrapolated by a common approach, based on 
the dark matter halo formation rate, calibrating the star-formation efficiency to the SFR at 
redshift z < 6. From this extrapolated SFR density, we can then obtain the intrinsic supernova 
rate by using the properties of core-collapse SNe in the local universe as a guide. 
In particular, we can estimate the global SFR density at redshift z as 



where e.^ is the efficiency (by mass) for the conversion of gas into stars, MdM{dn{> 
M, z) /dM) is the mass density contributed by halos with total (dark matter -I- baryonic) masses 
between M and M + dM, t{z) is the age of the universe at redshift z, and P{z) is the proba- 
bility per unit time that new stars form in a mass element of age r = t{z} —t{z') (normalized 
to dxP{t) = 1). We adopt the fiducial value of e» = 0.1 (see, e.g., EH)- Note that the star 
formation efficiency in minihalos (i.e. halos with virial temperatures below 10^ K) that contain 
pristine metal-free gas could be significantly lower than 0.1, as suggested by numerical simu- 
lations of the first generation of stars mi20l . The expected pre-reionization SNe rates would 
then lie closer to our T^^^ > lO'* K curves prior to reionization (see discussion below), making 
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the detection of the reionization feature significantly more difficult. However, the efficiency is 
likely to be very sensitive to even trace amounts of metallicity [21], and conditions for star- 
formation may result in a standard initial mass function (IMF) in gas that has been enriched 
to metallicities above a fraction 10^"* of the solar value. Indeed, it is unlikely that metal-free 
star-formation in minihalos can produce enough ionizing photons to dominate the ionizing 
background at reionization (e.g. 1721 1741 [157 1). We assume further that star-formation occurs 
on an extended time-scale, corresponding to the dynamical time, f^y^ ~ [Gp(z)]^'/^ I30II66I : 



P{'c) = -^exp 

'dyn 



T 

t, 



dyn 



(2) 



where p{z) ~ ZicPci-it(z) is the mean mass density interior to collapsed spherical halos (e.g. 
(9J), and is obtained from the fitting formula in L25J , with Ac = ISn^ ~ 178 in the Einstein- 
de Sitter model. The minimum mass, ^im^iz) in ec[. (1), depends on the efficiency with which 
gas can cool and collapse into a dark matter halo. Prior to reionization and without molec- 
ular hydrogen, corresponds to a halo with virial temperature, Tvii- ~ 10'* K; with a 
significant H2 abundance, the threshold decreases to Tyir 300 K (|72|; we use the conver- 
sion between halo mass and virial temperature as given in |9|). Post reionization, the Jeans 
mass is raised, so M^[„{z) could increase. The degree of self-shielding, the ability of the halo 
gas to cool, as well as the amount of H2 present in the high-redshift low-mass halos is un- 
certain 1361 , and so below we present results for several values of Minin(z), which we will 
henceforth express in terms of Ty[^(z). 

Next, from p*(z) we obtain the intrinsic differential SNR (number of core collapse SNe 
per unit redshift per year) with 

dN 1 dV(z) , , 

-r = '7SNT- ^P*z, (3) 

dz 1 + Z dz 

where the factor 1/(1 + z) accounts for time dilation, dV{z)/dz is the comoving volume in our 
past light cone per unit redshift, and r7sN is the number of SNe per solar mass in stars. For a 
fiducial Salpeter initial mass function (IMF), we obtain 77sn ~ 1 / 1 80 Mg ' , assuming all stars 
with masses 9 Mg) < M < 40 become core collapse SNe f841. We neglect the lifetime of 
these high mass stars in determining our SNR; this is a reasonable assumption as the lifetimes 
(as well as the spread in the lifetimes) are shorter than a unit redshift interval for redshifts of 
interest. Note that an alternative extreme shape for the IMF, consisting entirely of 100-200 
Mq stars fTll20l would yield a similar value for t)sn- 

We present our SFR densities {top panel) and SNRs {bottom panel) in Figure[T] The curves 
correspond to redshift-independent virial temperature cutoffs of T^{^ = 300, 10"*, 4.5 x 10^, 
and 1.1 X 10^ K (or circular velocities of Vcijc = 3, 17, 35, and 55 km s^', respectively), top 
to bottom, spanning the expected range I173II36I . Also shown are results from GOODS |64|: 
the bottom points assume no dust correction and the top points are dust corrected according to 
(3 1 ; the statistical error bars lie within the points fSl . As there are large uncertainties associated 
with dust correction, each pair of points (top and bottom) serves to encompass the expected 
SFR densities. 

Our SFRs are consistent with other theoretical predictions (e.g. |158| [8l l23l ). as well as 
other estimates from the Hubble Ultra Deep Field |26 | and the FORS Deep Field on the VLT 
(601, after they incorporate a factor of 5-10 increase in the SFR |3| due to dust obscuration. 
Furthermore, we note that our SNRs, which at z > 5 yield 0.3 - 2 SNe per square arcminute 
per year, are in good agreement with the ~ 1 SN per square arcminute per year estimated 
by 11221 by requiring that high-redshift SNe produce a mean metallicity of ~ 0.01 Zr.j by 
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z ^ 5. Note, however, that the rates we obtain are significantly higher (by a factor of ~ 60 
- 2000 at z ~ 20) than those recently found by 1 185 |. The reason for this large difference is 
that Wise & Abel consider star-formation only in minihalos, and they assume a very low star- 
formation efficiency of a single star per minihalo, as may be appropriate for star-formation out 
of pristine (metal-free) gas in the first generation of minihalos 1 1 '201. In contrast, we assume 
an efficiency of =0.1, which may be more appropriate for star-formation in pre-enriched 
gas that dominates the SFR just prior to reionization (including star-formation in minihalos). 

As mentioned above, by increasing the cosmological Jeans mass, reionization is expected 
to cause a drop in the SFR (and hence the SNR), with the rates going from the horizontally 
striped region in Figure[T]at z > Zre to the vertically striped region at z < Zre- 

The redshift width of this transition is set by a combination of large-scale cosmic variance, 
radiative transfer, and feedback effects. For the majority of the paper, we use Azie ~ 1 as 
a rough indicator of the width of the transition we are analyzing. We distinguish between 
"reionization" and a "reionization feature", and use Azre as an indicator of the width of the 
later. Even with an extended reionization history {Az ~ 10), fairly sharp (Azre < 3) features 
are likely, as discussed in detail in [1 19| . 

The other important factor determining the usefulness of the method proposed here is the 
factor by which the SFR drops during the reionization epoch. The size of this drop is medi- 
ated by the effectiveness of self-shielding and gas cooling during photo-heating feedback: i.e. 
on whether or not the star-formation efficiency is significantly suppressed in those halos that 
dominate the SFR and SNR immediately preceding the reionization epoch. Given the uncer- 
tainties about this feedback discussed above, we will consider a range of possibilities below, 
parameterized by the modulation in the virial temperature threshold for star-formation during 
reionization. 



6.3 The Rate of High-Redshift SNe Detectable in a Future Survey 

Given the intrinsic star-formation and SN rates, our next task is to estimate the number of SNe 
that could actually be detected in a future search with JWST. In general, the number of SNe 
per unit redshift, dN^xp/dz, that are bright enough to be detectable in an exposure of duration 
fexp can be expressed as 

dNexp dN f" 

—, — = ^ /sn(> tohs) dtohs , (4) 
dz dz Jo 

where {dN / dz)dta^^ is the number of SNe which occurred between r^bs and fob, +ii/obs ago 
(per unit redshift; note that the global mean SNR will evolve only on the Hubble expansion 
time-scale, and can be considered constant over several years), and /sn(> 'obs) is the fraction 
of SNe which remains visible for at least f^bs in the observed frame. Then the total number of 
SNe detected in a survey of duration rsurv is 

_ gsurv AQy'OV 
-''surv ~ ^ . A J'exp , yj) 

^ ^exp 

where ASl^aw is the instrument's field of view, and /sui-v/(2fexp) is the number of fields which 
can be tiled in the survey time, fsuiv (we add a factor of 1/2 to allow for a second pair of filters to 
aid in the photometric redshift determination; note that this provides for imaging in 4 different 
yVKST bands; see discussion below). Note also that equation ^ is somewhat idealized, in that 
it assumes continuous integration for a year, and e.g., does not account for time required to 
slew the instrument to observe different fields. In principle, each field has to have repeated 
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Fig. 1. Upper Panel: SFR densities obtained in a model based on dark matter halo abundances. 
The curves {top to bottom) correspond to different lower cutoffs on the virial temperatures of 
star-forming halos, Tyir > 300, 10'*, 4.5 x 10"*, and 1.1 x 10^ K (corresponding to circular 
velocity thresholds of v^ijc > 3, 17, 35, and 55 km s^'). Dots indicate results from GOODS 
1641 : the lower set of points assume no dust correction, while the upper set of points are dust 
corrected; statistical 1-ct error bars lie within the points. The lines connecting each pair of 
points span the expected range of SFR densities. Lower Panel: the total global supernova rate 
accompanying the SFR densities in the top panel (adapted from H 19i ). 

observations (to detect SNe by their variability), and therefore any dedicated survey should 
target fields that have already been observed. Furthermore, a dedicated year-long program 
may not be necessary, because the effect could be detected with relatively few fields (at least 
under optimistic assumptions; see discussion below), and several fields with repeated imaging 
(separated by > 1 - 2 years) may already be available from other projects; these fields can 
then be used for the SN search. 

In general, /sn(> 'obs) in equation l|4j, i.e. the fraction of SNe which remain visible for 
at least fgbs' depends on (i) the properties of the SN, in particular their peak magnitude and 
lightcurve, and the distribution of these properties among SNe, and (ii) on the properties of the 
telescope, such as sensitivity, spectral coverage, and field of view. In the next two subsections, 
we discuss our assumptions and modeling of both of these in turn. 
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Empirical Calibration of SN Properties 

At each redshift, we run Monte-Carlo simulations to determine /sn(> ?obs) in equation We 
use the observed properties of local core-collapse SNe (CCSNe) in estimating /sn(> fobs)- 
For the high redshifts of interest here, we only consider core collapse SNe of Type II. SNe 
resulting from the collapse of Chandrasekhar-mass white dwarfs (Type la) are expected to be 
extremely rare at high redshifts (z > 6), as the delay between the formation of the progenitor 
and the SN event (> 1 Gyr; 11671 ') is longer than the age of the universe at these redshifts. 
Local core CCSNe come in two important varieties, types IIP and IIL, differentiated by their 
lightcurve shapes. We ignore the extremely rare additional CCSN types, e.g Type Iln and lib, 
which appear to have significant interaction with circumstellar material and constitute less 
than 10% of all CCSNe. Type Ib/c, which may or may not also result from core collapse, have 
luminosities and light-curves that are similar to Type IIL and occur less frequently. While the 
relative numbers of Type IIP and Type IIL SNe are not known even for nearby SNe, estimates 
imply that they are approximately equal in frequency |27|. We therefore assume that 50% of 
the high-redshift SNe are Type IIP and 50% are Type IIL. 

CCSNe result from the collapse of the degenerate cores of high-mass stars. The luminosity 
of CCSNe is derived from the initial shock caused by the core-collapse which ionizes mate- 
rial and fuses unstable metal isotopes (see |102| and references therein for a more detailed 
description of SN lightcurves). In the early stages of the SN, the shock caused by the core col- 
lapse breaks out from the surface of the progenitor (typically high mass red giants), resulting 
in a bright initial peak in the light curve that lasts less than a few days in the rest-frame of the 
SN. As the shock front cools, the SN dims. However, the SN may then reach a plateau of con- 
stant luminosity in the light curve, believed to be caused by a wave of recombining material 
(ionized in the shock) receding through the envelope. The duration and strength of this plateau 
depends on the depth and mass of the progenitor envelope, as well as the explosion energy, 
with those SNe exhibiting a strong plateau classified as Type IIP. A typical plateau duration 
is < 100 rest-frame days 1 135 1. In Type IIL SNe, this plateau is nearly non-existent, and the 
lightcurve smoothly transitions from the rapid decline of the cooling shock to a slower decline 
where the luminosity is powered by the radioactive decay of metals in the SN nebula. After 
the plateau. Type IIP SNe also enter this slowly declining 'nebular' phase. 

These observationally determined behaviors have been summarized in a useful form as 
lightcurve templates in |40|. We use these template lightcurves in determining /sn(> 'obs)> 
and normalize the lightcurves using Gaussian-distributed peak magnitudes (i.e. log-normally 
distributed in peak flux) determined by 11401 from a large sample of local Type IIP and Type 
IIL SNe (see Table|2]l- We perform the Monte-Carlo simulations with both the dust corrected, 
and dust uncorrected values in 11401 , since the dust production history of the early universe is 
poorly understood and is essentially unconstrained empirically. 

We use a combined high-resolution HST STIS + ground-based spectrum of the Type IIP 
supernova SN1999em (the 'November 5th' spectrum of 1121 ) as the template SN spectrum in 
order to obtain A^-corrections (with the | 40| lightcurves given in the restframe B filter). This 
spectrum was obtained within 10 days of maximum light, i.e. during the initial decline of 
the SN brightness, and has been dereddened by Ay = 0.3 mag 1 12, 83 1. While the spectrum, 
and hence the A'-corrections, of SNe evolve during the lightcurve, this effect is not strong 
for the wavelengths of interest 11351 , especially since the lightcurve template we use is well 
matched to the wavelengths being probed by the observations we consider below (leading to 
small A'-corrections). We have also used this Type IIP SN template spectrum to calculate K- 
corrections for Type IIL SNe. This is necessary due to the lack of restframe UV spectra of Type 
IIL SNe that can be combined with optical spectra, and justifiable because the A^-corrections 
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are relatively small and the broadband colors of both Type IIP and IIL SNe (a measure of 
the spectral shape that determines the A'-corrections) are similar, at least in the optical II35I . 
Note that we assume that very-high redshift core-collapse SNe are similar to local SNe in 
their spectra, peak luminosities, and temporal evolution. However, these assumption do not 
significantly impact our conclusions below, as long as (i) the average properties of the SNe 
do not change rapidly at high redshift (which could mimic the reionization drop), (ii) they 
do not become preferentially under-luminous (which would make high-z SNe less detectable, 
lowering the statistical confidence at which the reionization drop is measured), and (iii) the 
detection efficiency does not change rapidly with redshift (e.g. due to instrument parameters 
or spectral lines). 

Table 2. Means and standard deviations of the adopted peak absolute magnitudes of core 
collapse SNe. Values are taken from II140I . Note: an Mb = — 17 SN would be detectable out 
to z PS 8.2 at the flux threshold of 3 njy in the 4.5 ^m band with JWST. 



Corrected for Dust Not Corrected for Dust 



SN Type 


(Mb) 


CT 


(Mb) 


tr 


IIP 


-17.00 


1. 12 


-16.61 


1.23 


IIL 


-18.03 


0.9 


-17.80 


0.88 



In order to use the method outlined above to probe reionization, the SN redshifts must 
also be known to an accuracy of Az< I. SN redshifts can be determined via spectroscopy of 
either the SN itself, or of the host galaxy. However, as we have already noted, the host galaxies 
may only be marginally detectable even in imaging, and the SNe may be too faint for anything 
other than extremely low resolution spectroscopy. We present here only a very brief example 
of the possibility of obtaining redshifts from the extremely low resolution (X/AX ~ 5) spectra 
provided by multi-band imaging: a complete investigation of this possibility is warranted, but 
is beyond the scope of this paper. 

To the extent that Type II SNe spectra can be represented as a sequence of blackbodies of 
different temperatures (e.g. 1341 ) photometric redshifts will be impossible to obtain without 
information about the SN epoch, since temperature and redshift would be degenerate. How- 
ever, local Type IIP SNe show significant deviations from a blackbody in the UV (A < 3500A) 
due to metal-line blanketing in the SN photosphere, providing spectral signatures that could 
be used as redshift indicators, depending on their strength. In Figure |2] we show the evolution 
with redshift of the infrared colors (in bands accessible with JWST; see below) of our template 
SN spectrum, compared with the color evolution of a blackbody. The figure shows that the 
template spectrum deviates significantly from a blackbody. If the spectrum of the SN is always 
the same as the template spectrum, then there are good prospects for obtaining photometric 
redshifts for these SNe, at least in the redshift range z = 7 — 13. While the figure shows that 
multiple redshifts may be possible at fixed observed colors, the degenerate solutions would 
correspond to z > 16 SNe; contamination from such high redshift will be mitigated by the fact 
that these SNe are likely to be too faint to be detected. Of course, more detailed studies of the 
UV behavior of local SNe, especially their variety and spectral evolution, will be necessary to 
confirm the possible use of photometric redshifts. 
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Fig. 2. Infrared colors of a Type II SN as a function of its redshift. Solid curve: The template 
spectrum of the text, from the Type IIP SN 1999em. Diamonds are placed at intervals of Az = 
1 , labeled with the redshift. Error bars of 0. 1 5 mag are shown to give a sense of the photometric 
errors that may be expected, excluding the many possible systematic effects. Dashed curx'e: 
The colors of a blackbody at different redshifts (temperatures) are shown for reference. The 
asterisk marks the color of a 9000K blackbody at z = 10, triangles are at intervals of Az= 1 ; or 
equivalently, at a constant redshift, but with correspondingly different temperatures (adapted 
from fTT9l ). 

SNe Detectability and Survey Parameters 

As a specific example for the number of detectable SNe in a future sample, we consider ob- 
servations by JWST. The relevant instrument on JWST is NIRcam0 a near-infrared imaging 
detector with a FOV of 2.3' x 4.6'. A field can be observed in two filters simultaneously. NIR- 
cam will have five broadband filters (with resolution X/AX^ 5). We model the filter response 
as tophat functions with central wavelengths of 1.5, 2.0, 2.7, 3.5, and 4.5 /xm. For concrete- 
ness, below we will present results only for the 4.5 and 3.5 /im filters, since they are the 
longest-wavelength iVVSr bands; however, we allow time for imaging in two other bands, if 
needed for photometric redshift determinations. The current estimate of the JWST detection 
threshold at 4.5 jlm is > 3 nJy for a 10 (J detection and an exposure time of 10^ s. The 3.5 
/xm band is more sensitive, with a detection threshold of > 1 nJy for a 10 tr detection and an 
exposure time of 10^ s. 

We show our results for /sn(> fobs) in Figure [5] The solid curves correspond to z = 7, 
the dashed curves correspond to z = 10, and the dotted curves correspond to z = 13. In each 
panel, the top set of curves assumes a flux threshold of 3 nJy (or an exposure time of ?exp — 

10^ 

seconds in the 4.5 /im band), and the bottom set assumes 9.5 nJy (fexp = 10** seconds in the 
4.5 jUm band, background dominated). The top panel further assumes no dust extinction, while 
in the bottom panel, we adopt the same dust extinction as in the low redshift sample f 1401. 
Understandably, the distributions get wider as redshift increases (due to time dilation), but 
the total visible fraction, /sn(> 0), gets smaller (due to the increase in luminosity distance). 



^ See |http ://ircamera. as . arizona. edu/ntrca m for further details. 
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Fig. 3. The fraction of SNe which remain visible for an observed duration of /^bs or longer. 
The curves correspond to SN redshifts of z = 7 (solid curx'e), z = 10 {dashed ciirx'e), and 
z = 13 {dotted curve). In each panel, the top set of curves assumes a flux threshold of 3 njy 
(exposure time of ?exp = 10^ s with the 4.5 /xm7W5rband), and the bottom set assumes 9.5 nJy 
(fexp = 10"* s). The top panel assumes no dust extinction, while the bottom panel assumes the 
same dust extinction as observed within the low-redshift sample 11401 (adapted from 11 191 ). 

The double bump feature in some of the curves corresponds to the plateau of Type IIP SNe 
lightcurves discussed above. 

As seen in Figure [3] most of the supernovae that are bright enough to be visible at all, 
will remain visible for up to --^ 1 - 2 years. Hence, in order to catch the most SNe, it will be 
necessary to have repeat observations of the SN survey fields a few years apart, to insure that 
most of the observed SNe will be identified as new sources or sources that have disappeared. 
A time between observations that is comparable to or larger than the SN duration will be the 
optimal strategy for detecting the most SNe 11261 . Note that we ignore the time required to 
collect reference images, since observations conducted for other programs will likely provide 
a sufficient set of such reference images. However, as shown in equation (O, we allow time 
for the field to be imaged in four different bands, to aid in photometric redshift determination. 

To be more explicit, we find that in order to obtain the largest number of high-redshift 
SNe, in general it is a more efficient use of JWSTf, time to 'tile' multiple fields rather than 
'stare' for extended periods of time (> 10^ s) at the same field f 1261. This is because of 
strong time dilation at these redshifts. As can be seen from Figure [3] most detectable SNe 
will remain above the detection threshold for several months, even assuming 10'* s integration 
times. Also evident from Figure[3]is that the increase in the total visible fraction of SNe going 
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from fexp = 10^ s to f^xp = 10^ s, is less than the factor of 10 increase in exposure time. As a 
result, a fiducial 1-yr JWST survey would therefore detect more SNe using Jgxp = 10"* s than 
using fexp = 10^ s (see Figures |4] and |5j- Understandably, this conclusion does not hold for 
very high-redshifts, z > 14, where SNe are extremely faint, and require very long exposure 
times to be detectable. However, even with such long exposure times, very few SNe will be 
detectable at these large redshifts, rendering the use of longer exposure times unnecessary. 

SNe Detection Rates 

The number of SNe that could be detectable in putative future surveys are shown in Figures 
|4]and|5] The curves correspond to the same virial temperature cutoffs for star-forming halos 
as in Figure [T] Solid lines assume no dust obscuration; dashed lines include a correction for 
dust obscuration as discussed above. Figure|4]shows results assuming flux density thresholds 
of 9.5 njy (or fexp = 10'* s with the 4.5 flm 71^57 filter) (top panel) and 3 njy (/exp = 10^ s) 
(bottom panel). Figure [5] shows results with the 3.5 flm filter assuming equivalent exposure 
times: flux density thresholds of 3.2 nJy (fexp = 10'* s with the 3.5 /xm yVKST filter) (top panel) 
and 1 nJy (fexp = 10^ s) (bottom panel). The right vertical axis displays the number of SNe 
per unit redshift per FOV (2.3' x 4.6'); the left vertical axis shows the number of SNe per unit 
redshift in a fiducial 1-year survey. As mentioned above, reionization should be marked by a 
transition from the region bounded by the top two solid curves to the region enclosed by the 
bottom three solid curves (or the analog with the dashed curves if dust is present at the time 
of reionization). 

We note that our expected rates are somewhat higher than those in 1341 . a previous study 
which included SNe lightcurves and spectra in the analysis. For example, we find 4-24 SNe 
per field at z > 5 in the 4.5 /im filter with fexp ~ 10^ s, compared to ~ 0.7 SNe per field at z > 5 
obtained by |34| (after updating their JWST specifications to the current version). However, 
they use SFRs extrapolated from the low-redshift data available at the time, which are not a 
good fit to recent high-z SFR estimates I64ll60|[26l , and are lower than our z > 5 SFRs by a 
factor of 6 - 400 Taking this factor into account, their procedure yields 4-27 SNe per field 
at z > 5, which is in excellent agreement with our estimate of 4 - 24 SNe per field at z > 5. 

Pop-Ill SNe 

As discussed above, high-redshift SNe whose progenitor stars are formed from metal-free gas 
within minihalos could be intrinsically very different from the low-redshift SNe, due to differ- 
ences in the progenitor environments (e.g. very low metallicities; 1 1 20 1). If such differences 
could be identified and detected, then these "pop-Ill" SNe could provide valuable information 
about primordial stars and their environments. Indeed, 1 185 1 studied the redshift-distribution 
of such primordial SNe, but only briefly addressed the issue of their detectability. In order 
to directly assess the number of such SNe among the hypothetical SNe samples we obtained 
here, in Figure [6] we plot the fraction of SNe whose progenitor stars are located in miniha- 
los. The solid curve assumes our fiducial model with a Salpeter IMF and £*rninihalo ~ 

0.1; the 

dotted curve assumes that each minihalo produces only a single star, and hence a single SN, 
over a dynamical time (assuming that strong feedback from this star disrupts any future star 

^ The possibility that the SFR in a "Lilly-Madau diagram" remains flat, or even increases, at 
redshifts z > 5, owing to star-formation in early, low-mass halos, is also expected theoret- 
ically (see, e.g., Fig.l in (23 1 and associated discussion). 
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Fig. 4. The number of high-redshift SNe detectable with the 4.5 flm JWST filter. The curves 
correspond to the same virial temperature cutoffs as in Figure[T] Solid curves assume no dust 
obscuration; dashed curves adopt dust obscuration in the same amount as observed in the low 
redshift SNe sample. The figure shows results assuming flux density thresholds of 9.5 njy (or 
'exp =10^*5 with JWST) {top panel) and 3 nJy (or ?exp = 10^ s with JWST) {bottom panel). The 
right vertical axis displays the number of SNe per unit redshift per field; the left vertical axis 
shows the number of SNe per unit redshift found in ?surv / (2^exp) such fields (i.e. the differential 
version of eq. l|5j with fsurv = 1 yr)- Reionization should be marked by a transition from the 
region bounded by the top two solid curves to the region enclosed by the bottom three solid 
curves (or the analog with the dashed curves if dust is present at the time of reionization). 
Adapted from | fTT9l . 

formation; as in |185))- The figure shows that with an unevolving star formation efficiency, 
progenitor stars in minihalos would account for > half of the SNe at z > 9. On the other hand, 
in the extreme case of a single SNe per minihalo, progenitor stars in minihalos would account 
for less than ~ 1% of the z ~ 10 SNe (although at the earliest epochs, z > 22, they would still 
constitute > half of all SNe). Given the overall detection rate of several hundred z > 10 SNe 
with a 1 year JWST survey we have found above, even in this extreme case. Figure |6] implies 
that several of these SNe could be caused by pop-Ill stars. 
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Fig. 5. Same as Figure ID but with the 3.5 flm filter, instead of the 4.5 fxm one. Note that 
we present results for comparable exposure times, hence the sensitivity thresholds are three 
times lower than in Figure|4l due to the disparate sensitivities of the 3.5 ^m and 4.5 fxm filters 
(adapted from L119J ). 

6.4 Detecting Features due to Feedback 

In summary, we find that 4-24 SNe may be detectable from z > 5 at the sensitivity of 3 
njy (requiring 10^ s exposures in a 4.5 /xm band) in each ~ 10 arcmin^ JWST field. In a 
hypothetical one year survey, we expect to detect up to thousands of SNe per unit redshift 
at z --^ 6. These rates are high, and, if reionization produces a fairly sharp features in the 
reionization history (with a drop in the SFR by a factor of a few, spread at most over Azre ~ 
1 - 3), we have also shown 1 119| that the number of SNe is sufficient to detect the feature 
out to z ~ 13, as well as set constraints on the photo-ionization heating feedback on low- 
mass halos at the reionization epoch. Specifically, for a wide range of scenarios at Zre ^13, 
the drop in the SNR due to reionization can be detected at S/N > 3 with only tens of deep 
JWST exposures. These results therefore suggest that future searches for high-z SNe could be 
a valuable new tool, complementing other techniques, to study the process of reionization, as 
well as the feedback mechanism that regulates it. 
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Fig. 6. The fraction of SNe whose progenitor stars are located in minihalos (with virial temper- 
atures of Tvii- < lO'^K). These SNe may be pair instability SNe, which would be much brighter 
than the normal core-collapse SNe used in the computation of SN rates in the previous figures. 
The solid curve assumes our fiducial model with a Salpeter IMF and £*niinihalo — 0.1; the dot- 
ted curve assumes that each minihalo produces only a single star (and therefore a single SN) 
over a dynamical time (adapted from 1 1 19j). 



7 The Dark Ages as a Probe of The Small Scale Power Spectrum 

Since, as discussed above, GRBs are a rare tracer of the cosmic SFR, a mere presence of a 
GRB at, say, z > 10 will indicate that non-linear structures exist already at this redshift: the 
stars that give birth to the GRBs must form out of gas that collected inside dense dark mat- 
ter potential wells. Structure formation in a cold dark matter (CDM) dominated universe is 
"bottom-up", with low-mass halos condensing first. In the current "best-fit" cosmology, with 
densities in CDM and "dark energy" of {£2yi,£2/^) ~ (0.3, 0.7) that has emerged from WMAP 
and other recent experiments, DM halos with the masses of globular clusters, lO^^^Mp, con- 
dense from ~ 3a peaks of the initial primordial density field as early as z --^ 25. It is natural 
to identify these condensations as the sites where the first astrophysical objects, including 
the first massive stars, were born. As a result, one expects to possibly find GRBs out to this 
limiting redshift, but not beyond. 

With a scale-invariant initial fluctuation power spectrum, the CDM theory has been re- 
markably successful, and matched many observed properties of large-scale structures in the 
universe, and of the cosmic microwave background radiation. However, the power spectrum 
on scales corresponding to masses ofM < 10'' Mp, remains poorly tested. Some observations 
have suggested that the standard model predicts too much power on small scales: it predicts 
steep cusps at the centers of dark matter halos, whereas the rotation curves of dwarf galaxies 
suggest a flat core; it also predicts more small satellites than appear to be present in the Local 
Group (reviewed by, e.g. f 1491). Although astrophysical explanations of these observations 
are possible, several proposals to solve the problem have been put forward that involve the 
properties of dark matter. These include self-interacting dark matter |T61 1, adding a repulsive 
interaction to gravity I69II136I , the quantum-mechanical wave properties of ultra-light dark 
matter particles 1871 . and a resurrection of warm dark matter (WDM) models 1 18J . 
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By design, a common feature of models that attempt to solve the apparent small-scale 
problems of CDM is the reduction of fluctuation power on small scales. The loss of small- 
scale power modifies structure formation most severely at the highest redshifts, where the 
number of self-gravitating objects is reduced. In each model, there exists a redshift beyond 
which the number ofGRBs (or any other object) is exponentially suppressed; a detection of a 
GRB beyond this redshift can be used to rule out such models. As an example, ref.|6| showed 
that in the case of WDM models, invoking a WDM particle mass of ~ IkeV (the mass needed 
to solve the problems mentioned above), it becomes difficult to account for the reionization of 
the universe by redshift z ~ 6, due to the paucity of ionizing sources beyond z = 6. 

GRBs, if discovered in sufficient numbers at z > 6, have the potential to improve this 
constraints significantly, since increasingly higher redshifts probe increasingly smaller scales. 
The constraints, in particular, that the detection of distant GRBs would place on structure 
formation models with reduced small-scale power was quantified in [ 120|. In this work, we 
computed the number of GRBs that could be detectable by the Swift satellite at high redshifts 
(z > 6), assuming that the GRBs trace the cosmic star formation history, which itself traces 
the formation of non-linear structures. Simple models of the intrinsic luminosity function 
of the bursts were calibrated to the number and flux distribution of GRBs observed by the 
Burst And Transient Source Experiment (BATSE). Under these assumptions, the discovery of 
high-z GRBs would imply strong constraints on models with reduced small-scale power. For 
example, a single GRB at z > 10, or, alternatively, 10 GRBs at z > 5, discovered within a 
two-year period, would rule out an exponential suppression of the power spectrum on scales 
below = 0.09 Mpc (exemplified by warm dark matter models with a particle mass of m^ = 2 
keV). Models with a less sharp suppression of small-scale power, such as those with a red tilt 
or a running scalar index, n.,, are more difficult to constrain, because they are more degenerate 
with an increase in the power spectrum normalization, cTg, and with models in which star- 
formation is allowed in low-mass minihalos. We find that a tilt of Sus ~ 0. 1 is difficult to 
detect; however, an observed rate of 1 GRB/yr at z > 12 would yield an upper limit on the 
running of the spectral index, a = dns/dlnk > —0.05. 
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